Long-range and frustrated spin-spin interactions in crystals of cold polar molecules 
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We describe a simple scheme for the implementation and control of effective spin-spin interactions in self- 
assembled crystals of cold polar molecules. In our scheme spin states are encoded in two long-lived rotational 
states of the molecules and coupled via state dependent dipole-dipole forces to the lattice vibrations. We show 
that by choosing an appropriate time dependent modulation of the induced dipole moments the resulting phonon- 
mediated interactions compete with the direct dipole-dipole coupling and lead to long-range and tunable spin- 
spin interaction patterns. We illustrate how this technique can be used for the generation of multi-particle 
entangled spin states and the implementation of spin models with longe-range and frustrated interactions which 
exhibit non-trivial phases of magnetic ordering. 

PACS numbers: 03.67.Lx, 37.10.Pq, 75.10.Jm 



I. INTRODUCTION 

The quantum Ising model describes a system of interact- 
ing spins where the coupling among the spins competes with 
a transverse magnetic field. It is one of the simplest models 
which captures many essential aspects of quantum magnetism 
and has been successfully applied to study the transition from 
a paramagnetic phase to ferro- or anti-ferromagnetic ordering 
01 • Due to its practical relevance this model has been widely 
studied and for many realizations with nearest neighbor inter- 
actions the ground and thermal states predicted by this model 
are well understood. However, less is known about more gen- 
eral Ising models with long-range or frustrated spin-spin inter- 
actions [2], where already identifying the ground state prop- 
erties can be a complicated and numerically demanding task. 
Apart from its important role in condensed matter physics the 
Ising and related spin models have recently also attracted a 
lot of attention in the field of quantum information process- 
ing. Here the dynamics generated by a controlled Ising inter- 
actions between qubits can be used to generate various two 
and multi-partite entangled states as a fundamental resource 
for quantum computation 0-0]. These potential applications 
have over the past years stimualted a lot of additional work 
in the field of spin models focused in particular on the out 
of equilibrium dynamics and entanglement properties of this 
system 

The broad interest in a more detailed understanding and 
control of interacting spin systems on a small and larger scale 
has stimulated various proposals for the implementation and 
simulation of quantum Ising models using isolated atomic or 
coherent solid state systems |7|]. A pioneering role in this 
context is played by systems of trapped ions [|8l |9[], where 
effective, phonon-mediated spin-spin interactions can be im- 
plemented and controlled by applying state-dependent light 
forces | llOt4l6ll . Indeed, several proof-of-principle experi- 
ments IU7I420I1 have already demonstrated the possibility to 
simulate Ising interactions with up to nine spins and with such 
ion trap quantum simulators it might soon be possible to out- 
perform the best numerical simulations which are achievable 
on classical computers today. Based on recent advances in 



cooling and trapping techniques for diatomic molecules l2~ll - 
l28l . it is expected that in the near future a similar level of 
control can be achieved with ensembles of ultra-cold polar 
molecules. In this system spin states can be encoded in long- 
lived rotational or hyperfine states of the molecules, which can 
then be manipulated with microwave fields and coupled via 
strong electric dipole-dipole interactions ll29j[30l . Compared 
to ions, polar molecules can easily be trapped in standard op- 
tical lattice potentials of different geometries or - by aligning 
their dipole moments - be stabilized in a high density crys- 
talline phase ll3Tl - [34l . The combination of these exceptional 
properties make polar molecules one of the most promising 
systems for larg e scale simulation of non-trivial quantum spin 
models flUS. 



In a recent work H43H we have analyzed the implementa- 
tion of spin-spin interactions in a system of polar molecules 
prepared in a self-assembled dipolar crystal under ID or 2D 
trapping conditions. We have shown that in this setting reso- 
nantly enhanced phonon-meditated spin-spin interactions can 
dominate over the direct couplings and we proposed to use 
this feature for the implementation of local qubit operations 
in molecular quantum computing schemes IH41 l45l . In this 
paper we extend this analysis and investigate potential appli- 
cations of this technique for the design of more general Ising 
models with non-trivial spin-spin interaction patterns. In par- 
ticular, we show that in close analogy to trapped ion system 
we can use a time-dependent modulation of the induced dipole 
moments to address collective phonon modes in the crystal. 
The competition between short-range direct and long-range 
phonon-mediated spin-spin interactions then allows us to tai- 
lor the resulting effective spin-spin couplings and tune the in- 
teraction strength and signs freely via changing the detuning 
between the dipole moment and phonon frequencies. We il- 
lustrate how this technique can be used in the case of a small 
ID crystal of polar molecules to implement frustrated spin- 
spin couplings and discuss applications for the generation of 
multi-particle entangled states. Further, we provide an out- 
look, how this technique could be extended for simulation of 
larger spin systems. 

The remainder of the paper is organized as follows. In Sec. 
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FIG. 1. (Color online) Setup, (a) Polar molecules with two internal 
states |g) and \e) are confined in a ID tube by a strong transverse 
trapping potential. An electric dc field £ dc is used to align the dipole 
moment of the polar molecules while an additional weaker ac field 
£ a c(0 is used to induce a time dependent modulation of dipole mo- 
ment, which is of opposite sign for the two states as shown in (b). 
See text for more details. 



II we introduce our model and present a general derivation of 
the resulting Ising Hamiltonian. In Sec. Ill we discuss the 
validity of our spin model and how it can be used to generate 
a graph state using a simple toy system of three molecules in a 
harmonic trap. In Sec. IV we use the same setting to illustrate 
the frustrated character of this Ising Hamiltonian and describe 
the accessible 'phases' and the entanglement properties of the 
ground states. In Sec. V we discuss potential generalization 
of this techniques to larger spin systems and finally in Sec. 
VI we present more details on the encoding of qubits states 
in rotational states of the molecules. The main results and 
conclusions of this work are summarized in Sec. VII. 



the molecules to encode our effective spin states and denote 
them by \g) and \e). The corresponding induced dipole mo- 
ments {g\fj\g) - fi g e z and (e\f}\e) = fi e e z are a function of the 
applied bias field Eb and we assume that the dipole moments 
satisfy the following conditions: i) fJ. e (Eo) = /-i g (Eo) = /j for 
a specific value of the bias field Eq, ii) Oe^Eo = — dgfig^ 
and iii) «: fio. In Sec. VI below we will describe in 

more detail how these conditions can be achieved for a spe- 
cific set of rotational states and related ideas have been dis- 
cussed in Ref. 01 IH @1. Under the validity of assump- 
tions i)-iii) the combination E/,(t) = Eq + E ac cos(wf) of a 
static and a weak oscillating bias field then allows us to imple- 
ment a dipole operator of the form jj = fioe z (l + £cos(ojt)cr ), 
where e <K 1 and cr = \g)(g\ - \e){e\ is the Pauli operator. 
We write the dipole-dipole interaction as a sum of two parts, 
V dd (xj - Xj) = V° dd + V l dd (t) where for D = fj.l/(4ne ) the static 
and state independent part is given by 



V 



while the remaining part 
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depends on the oscillating dipole moments and couples inter- 
nal and external degrees of freedom. 



II. SPIN-SPIN INTERACTIONS IN MOLECULAR 
DIPOLAR CRYSTALS 



We consider a set of N polar molecules with dipole mo- 
ments aligned along the z-axis by an external bias field Eb = 
Ete z and their motion confined to the x,y plane by a strong 
optical or electric trapping potential. For simplicity we will 
focus in the following mainly on a one-dimensional (ID) trap- 
ping configuration as shown in Fig. [1] but our results can be 
generalized to 2D or intermediate trapping conditions. For 
molecules of mass m the system is then described by the 
Hamiltonian 
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where x, and are the position and momentum operators of 
the molecules and v denotes the trapping frequency of an addi- 
tional weak confinement potential along the tube. The dipole- 
dipole interaction is 

„ n 1 V Mj ~ 3(M)(M) 

V dd ({xi}) = > ■ , (2) 

8tt6 j^j m-xf 

where the /?, is the dipole operator of the z'-th molecule. 

In the following we chose two long-lived internal states of 



A. Spin-phonon interactions in a dipolar crystal 



In the limit e <k 1 the dipole-dipole coupling is dominated 
by the state independent part V dd given in Eq. (f3]l which in 
analogy to the Coulomb interaction in trapped ion systems 
IE 0] stabilizes the molecules against close encounter col- 
lisions and leads to the formation of a dipolar (quasi) crys- 
tal at low temperatures l31r - [34ll . In this limit molecules be- 
come localized at equilibrium positions x? and up to sec- 
ond order in the remaining small fluctuations 6x, the crys- 
tal dynamics can be described by the phonon Hamiltonian 
H p = H(e — > 0) = 2«^ w « a n a « USS^]. Here al(a n ) and a>„ 
are the creation (annihilation) operators and phonon frequency 
of the nth mode respectively, and the total Hamiltonian (Q]) can 
now be written as 



H = H p + V dd (f). 



(5) 



In contrast to V ( ! , 

dd 



the term V l dd it) 



(en 2 , eVo) con- 
tains spin-dependent dipole-dipole interaction, but is reduced 
by e «; 1. Therefore, we can expand V\,(f) in terms of the 
small parameters e and \Sxi\/a, where a is the typical lattice 
spacing. To the lowest relevant order in these two parameters 
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we obtain |34|43|,|4 
~2 



V d/0 ~ I^^^J COs2 ( Wf ) + X n SnjCrj( a n + Cll) COS(Wf), 

(6) 



where 
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De 2 



i J' 



n ^ 3De(x° - x°) 
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(c„j - c„j). (8) 



Here the c„j are the normalized mode function amplitudes de- 
fined by SXi = Yjn ^K^ m( ^n)c„,i(a„ + a' n ) and in Eq. © we 
have already neglected an oscillating single spin term, which 
does not give a relevant contribution in the dynamics dis- 
cussed below. The first term in Eq. (O is the direct spin-spin 
interaction and the last term couples the internal states of the 
polar molecules with their external motion. In the interaction 
picture with respect to H p , the total Hamiltonian of the system 
simplifies to 

„ hG°. 

H(t) = 2j -^cr;o~COS(arf) 2 

i*i 

+ J] hg nJ cos(&*)<7f(a»e~ W + ale^-'). (9) 
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FIG. 2. (Color online) The normalized spin-spin coupling strength 
Gjj for three polar molecules in a ID harmonic trap as a function 
of the dipole moment oscillation frequency uj. The normal mode 
frequencies are indicated by the solid dots in the x-axis. 



HI. TOY SYSTEM : THREE MOLECULES IN A 
HARMONIC TRAP 

As a first example we will show in this section, how the 
general approach outlined above can be used to create a spe- 
cific spin entangled state, namely a so-called graph state, for 
the simplest case of three molecules in a harmonic trap. We 
use this example in particular also to discuss the achievable 
fidelities of such a state and the validity of our effective spin 
model. 



If we assume that the frequency of the dipole moment is not 
resonant with any phonon mode i.e., the condition |A„| = 
\a> - (jl>„\ » g„j, ~ii,n is satisfied, the excitations of real 
phonons can be avoided d ED, SMl IH gl On a 
timescale which is long compared to |A„| _I the phonon de- 
grees of freedom can be eliminated and the remaining slowly 
spin dynamics can be described by an effective Hamiltonian 
(see App.lAl 
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where G u = (G?. + G?.)/2 and 
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Note that in Eq. ( fTUb the bare spin-spin coupling G°. cos 2 (wf) 
has been averaged to G"/2. 



We see that in the effective Ising model in Eq. < TT~0b the 
phonon-mediated coupling Gjj is added to the bare dipole- 

dipole interaction G" which depending on the phonon struc- 
ture and the modulation frequency a> can in general lead to 
long-range and more complicated spin-spin interaction pat- 
terns. 



A. Harmonic trap 

We consider a set of N polar molecules which are trapped 
in a harmonic potential with trapping frequency v, such that 
the classical equilibrium positions x® are determined by the 
condition 



™^°-3Z)£ u(i 
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xi — x 



.0|5 



= 0, 



(12) 



and the distance between particles is in general not constant. 
We denote by a — min{|jc° - x°j\} the typical lattice spacing at 

the center of the trap which we write as a — ^ -\/D/{mv 2 ), 
where ij^ is a numerical constant which is £2 = ^6 for 
two molecules, £3 st 1.26 for three molecules and scales as 
£n ~ N~ 2 I 5 for larger N Q4I1 . By introducing dimensionless 
positions x, = Xj/a, the couplings normalized to the relevant 
dipole-dipole interaction energy can then be written as G° = 



fiG ./(D6 2 /« 3 ) 



gn,i 



SO 1 3 



~- 1/1*° 

a> 2 ) where a>„ - a>„/v and 
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(13) 



In Fig. [2] we plot the total spin-spin interaction strength G,y 
for N = 3 as a function of the dipole modulation frequency 
a). When we tune ai away from a specific phonon frequency 
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the phonon mediated spin-spin coupling contains contribu- 
tions from all modes except the center-of-mass (COM) mode 
where all c„j — 1 / V^V are the same. It shows that we can con- 
trol the coupling strength and sign completely by varying the 
detuning. One specific mode will dominate when we tune 
a) close to resonance. Since such an vibrational eigenmode is 
shared by all the particles this situation leads to long-range in- 
teractions and provides us with the opportunities to generate 
highly entangled spin states |6[]. In the opposite case of far 
detuning from all the normal modes, phonon-mediated spin- 
spin couplings will be strongly suppressed, and thereby the 
bare, almost nearest-neighbor interactions dominate. 



B. Generation of graph states 

In the following we illustrate how the effective spin-spin 
interaction entangles particles by presenting a numerical sim- 
ulation for the case of three trapped polar molecules. As 
shown in Fig. 12 for u> = 2.977v all the couplings between 
different particles are the same, i.e., G\2 — G23 — G13 — 
0.91 lDe 2 /(for 3 ). This corresponds to the Ising Hamiltonian 
H - TiijhGo^.o^j, which can directly be used for generation 
of a three-qubit graph state (D-Ht] 

l<5 «> = 23/2^1 + |e>l£^) ® (L?>2 + k> 2 cr|) ® (|g> 3 + |e> 3 crp, 

(14) 

by evolving an initial state (\g) + |e)) l8>3 /2 3 ^ 2 for a time t = 
tt/(4G 12 ). 

To see how the graph state can be generated and also to 
confirm the validity of our model, we compare the effective 
spin evolution of Eq. dTOb with the dynamics of the full time 
dependent Hamiltonian ©, which is evaluated in App. [A] 
Fig. 3(a) shows the evolution of the graph state fidelity 
P gmph (t) = Trp{p(f)|<D s ><<D s |} for e = 0.05 and e = 0.1 (solid 
lines), where p(f) is the density matrix evolved under the full 
Hamiltonian (O and Tr p denotes the partial trace over the 
phonon degrees of freedom. The dashed line shows the evolu- 
tion calculated from the effective spin Hamiltonian (TTOt . The 
plot clearly shows that our effective spin Hamiltonian agrees 
well with the full time dependent model and when Gnt — n/A 
three molecules are in the graph state with a fidelity of 96.5% 
in the case e = 0.1. For molecular dipole moments of a 
few Debye and lattice spacing of a few hundred nanometers 
the typical dipole-dipole interaction strength can be tens of 
kHz, in principle limited by the transverse optical trapping 
frequency only. If we set e — 0.1 we obtain G, ; /(2^) ~ 1 kHz 
and the typical time scale for the generation of the graph state 
is about 150 /is. Larger interactions could be achievable in 
electrostatic traps 15111 . where a stronger transverse confine- 
ment would allow dipole-dipole couplings D/a 3 in the MHz 
regime JM]. 

Discrepancies between the exact results and the effective 
Hamiltonian arise from higher order terms in the expansion 
parameters, fast oscillating terms which we have omitted in 
the effective spin dynamics and the possibility of phonons be- 
ing excited. The small high frequency oscillations superim- 
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FIG. 3. (Color online) (a) The occupation probability -P grap h of graph 
state dl4t calculated with the full Hamiltonian (dashed red line) and 
the effective spin Hamiltonian with e = 0.05 (solid green line) and 
0.1 (solid blue line), respectively, (b) The purity of the reduced state 
of the effective spin model as a function of e at time t = n/(4Gn)- 
In both plots we have assumed a> = 2.977v such that the nearest- 
neighbor interaction equal to the next-nearest-neighbor interaction, 

G\2 = G13. 



posed on the solid curve come from off-resonant couplings be- 
tween the spin and phonons and the fast oscillating bare cou- 
pling term cos{2a>t)o- z cr : . To study the parasitic entanglement 
between spins and phonons we plot in Fig. [2b) the purity of 
the reduced density matrix of the spin chain P s = Tr[p 2 ] for 
the final time t = 7r/(4Gn) where p s = Tr p [p] is the reduced 
spin density operator. The purity decays with e 2 , but specific 
values where all phonon modes rephase (maxima of P s ) can 
be exploited to achieve high fidelities for not too small cou- 
pling. In this plot we have assumed the the phonon modes 
are initially in the ground state. Actually, because the cou- 
pling between spins and phonons is off-resonant, the phonon- 
mediated spin-spin coupling is robust to thermal excitation of 
those modes [ 10], but for larger temperatures the fidelity will 
decay as ~ e 2 «,/, where «,/, m IcbT/Hv is the characteristic 
phonon occupation number. 

From our analysis we conclude that our effective Ising 
model given in Eq. ( fTUb is a valid approximation to the full 
system dynamics under the relevant conditions. This Ising 
model with long-range coupling can be used to generate 
highly entangled states, which are almost decoupled from the 
external motion of the molecules. Although cooling polar 
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region 


coupling strength 


ground state 


spin order 


transition 


la 
lb 


G12 > 0, G12 > G13 


lege), \geg) 


AFMS 


AFMS <-> FM 


Gi2,G 13 <0 


\eee), \ggg) 


FM 


II 


G12 > G13 > 


lege), \geg) 


AFMS 


AFMS <-> AFMA 


G13 >Gi2>0 


\eeg), \gee), \egg), \gge) 


AFMA 


III 


G12 > 0, Gn < 


lege), \geg) 


AFMS 





TABLE I. Summary of the spin order of the ground state of the three-spin system in the limit B* — > and for the regions I, II and III of the 
modulation frequency oj as described in the text. 



molecules is more difficult than ions, temperatures of T ~ 
mK and below have been reached in several recent experi- 
ments 023I428I1 . This temperature is sufficient to avoid un- 
wanted thermal effects on the effective spin-spin couplings 



IV. FRUSTRATED SPIN-SPIN INTERACTIONS 

In the last section we have discussed how to realize cr^cr-' 
interactions with different spatial profiles through exchanging 
virtual phonons. We now consider a simple extension of this 
model by adding a transverse field B x , 

H = %B X 2 of + J] fcGyof4 (15) 

i <<j 

This effective transverse field can be implemented by coupling 
the states \g) and \e) either directly via a resonant microwave 
field or via optical Raman transitions. For a homogenous 
system with nearest neighbor couplings the Ising model in 
Eq. ( TT3T > exhibits a well known transition from a paramagnetic 
to a ferromagnetic (FM) or antiferromagnetic (AFM) phase, 
depending on the sign of the G,,,+i [1]. For long-range spin- 
spin interactions the ground state of the Ising model is in gen- 
eral more involved as, for example, the competition between 
nearest and next-nearest neighbor coupling can lead to frus- 
tration effects where for B x — > the energy is minimized by a 
superpositions of multiple degenerate spin configurations |2|. 



A. Frustrated three-spin models with polar molecules 

Let us continue with the three spin system described above, 
which is the minimal setting where frustration effects can oc- 
cur 11911 . As mentioned above for B x — > and Gy < the 
energy is always minimized when all spins are aligned, i.e., 
\ggg) and \eee) and the ground state is ferromagnetic. For 
positive couplings the spins have to be anti-aligned to lower 
the energy and the ground state is AFM. It is possible to get 
the optimal configuration for a ID spin array which just has 
nearest neighboring interactions. But for a three-spin-array 
with long-range interactions, which can be seen as a 2D tri- 
angular lattice, it is impossible to arrange the spins in a way 



such that each spin is anti-parallel to all its interacting part- 
ners. The system then does not have a simple periodic ground 
state and is said to be frustrated [2]. For example, in our 
system when G12 or G13 > 0, the system is frustrated and 
the Ising ground state is an entangled superposition of six 
AFM states fl9tl : four antiferromagnetic asymmetric (AFMA) 
states \gge), \eeg), \gee), \egg) and two antiferromagnetic sym- 
metric (AFMS) states \geg), \ege) 0. However, the compe- 
tition between G12 and G13 will break the symmetry in the 
AFM order. For instance when G13 > G12 > 0, the next near- 
est neighbor interaction wins so that the spin ground state is 
the superposition of four degenerate AFMA states. 

As we discussed in Sec. Ill, the relative strengths and signs 
of the couplings Gh are controlled by the modulation fre- 
quency a> and as we can see from Fig. |2l this allows us to 
access several distinct parameter regimes depending on the re- 
lation between G12 and G13. In the following we denote by re- 
gion I, II and III the tree different regimes for to which are sep- 
arated by the two phonon mode frequencies 102 and W3. More 
precisely, for our model to be valid, a> can not be too close to 
resonance and has to satisfy the condition \a> - a>„\ » gn,;/2 
for any i and n. In the following we set e — 0.1 and re- 
strict the detunings to \a> - a> n \ > 10 x max{g n j/2}. Then 
the valid regions for to are: region I: < a) < 1.84v, region 
II: 2.63v < at < 3.23v and region III: co > 3. 78 v. In region 
I a change of the sign of both coupling constants occurs at 
u> = 1.35v and therefore we distinguish further between re- 
gion la and lb. Table U summarizes the ground states and spin 
order of the three molecules system in the five different pa- 
rameter regions in the limit B x — > 0. Since the parameters G,y 
change continuously in allowed regions, and the spin order 
changes at the points of G12 = (region I) and G12 = G13 (re- 
gion II), conventional 'phase transitions' happen even when 
B x = 0. When a> crosses a phonon mode (except the COM 
mode), Gjj change discontinuously so that the symmetry of 
the spins jumps as expected between different regions 111611 . 
However, these transitions are not accessible in this setting 
due to detuning requirements mentioned above. 

Fig. [4] summarize the 'phase diagram' of the three molecule 
system for a finite B x . To characterize the ground state we fol- 
low the approach used in Ref. [ 19] and plot in Fig. 4(a) the 
probability of both FM spin combinations P ggg + P eee . We see 
that when B x = the spin system is FM ordered in region lb 
and there is a phase transition between AFM order and FM 
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FIG. 4. (Color online) Ground state of the Ising model U5\ im- 
plemented with three molecules in a harmonic trap. The different 
regions I, II and III for the modulation frequency a> as defined in 
the text are separated by the phonon mode frequencies, (a) Phase 
diagram showing the overlap of the ground state with the two FM 
spin configurations, P ggg + P eee , as a function of B x /G ms where 

G lms = -\j(G 2 l2 + G| 3 + G| 3 )/3. (b) Phase diagram showing the 
overlap of the ground state with the two AFMS spin configuration, 
P + P 

1 eve • 1 pep- 



order in region I. To distinguish between AFMS and AFMA 
ground states we plot in Fig. 4(b) the population of two AFMS 
ground states, P ege + P geg - From this plot we can clearly iden- 
tify a transition between AFMS order and FM order in region 

I and a transition between AFMS and AFMA order in region 

II as indicated in Table U 



B. Adiabatic preparation of entangled ground states 
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FIG. 5. (Color online) Adiabatic preparation of the ground state 
of the three spin Ising model starting from an initial state |t/K0)} = 

| ) and reducing the transverse field according to B x (t)/G Tms = 

10 exp(-(vf) 2 /(507r)). The plots (a) and (b) show the overlap of \\ji{t)) 
with the states |AFMS) and |AFMA> for the values of to = 2.65v and 
a) = 3.2v, respectively. The plots (c) and (d) show the overlap of 
\i//(f)) with the entangled GHZ and W states as defined in the text. For 
these two plots the coupling parameters corresponding totii= 1.75v 
(c) and a> = 2.977v (d) have been assumed. 



To see the change in the spin order more directly, for in- 
stance the spin order changeing from AFMS to AFMA in re- 
gion II, we can look at the time evolution of the spin system 

when the initial state is | ) where |— ) = (\g) - \e))/ V2, 

which is the paramagnetic ground state for B x » G rms = 

^j(G 2 l2 + Gj 3 + G2 3 )/3. By adiabatically reducing B x /G rms 
the spin system will remain in the ground state of Eq. ( fl3l ). 
In Fig. 5 (a) and (b) we numerically evaluate the resulting 
time evolution for u> = 2.65v and a> = 3.2v respectively and 
plot the overlap of the evolved initial state with an equal su- 
perposition of two AFMS states, 

\AFMS) = Qgeg)-\ege))/y/2, (16) 

as well as the equal superposition of four AFMA states, 

| AFMA) = (\gge) + \egg) - \eeg) - \gee))/2. (17) 

We see that when B x /G rms 0, the population of the state 
| AFMS) is almost equal to 1 and the spin system shows AFMS 
order, while if the system in AFMA order, the final state would 
be | AFMA). 

While so far we have consider the adiabatic time evolution 
as an alternative way to identify the ground state phases of this 
system, such an experiment can also be of practical interest 
when the resulting ground state is an entangled state. From the 
discussion above we see that when we initialize each spin par- 



allel to a strong transverse field and then adiabatically lower 
the field compared to the Ising couplings, the adiabatic evo- 
lution of the initial state will result in an equal superposition 
of all Ising ground states and should therefore be entangled. 
For instance, in region lb, the spins show FM order and for 
a small but non-zero B x the actual ground state would be the 
GHZ state |GHZ) = (\ggg) - \eee))l v5. For the case of G n = 
Go > 0, the resulting ground state would be the equal super- 
position of six AFM states, which is a superposition of two W 
state, W = (\gge) + \egg) + \geg)-\eeg)-\gee)-\ege))l ^ H- 
The adiabatic preparation of these two special entangled states 
is summarized in Fig.|5]where in (c) we have set a> = 1 .75v to 
prepare the GHZ state and in (d) a> = 2.977v to evolve into a 
W superposition state. 



V. SCALABILITY 

In our discussion so far we have focused on a system of 
three molecule in a harmonic trap as the conceptually and 
experimentally simplest system where frustration phenomena 
can occur. In particular, a direct comparison with the anal- 
ogous trapped ions systems fl9tl shows that all the relevant 
parameter regimes of the frustrated Ising model are equally 
well accessible with polar molecules using a simple oscilla- 
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tion electric field. However, also in analogy to trapped ion sys- 
tems our techniques relies on a frequency selective addressing 
of individual phonon modes and is limited to molecular crys- 
tals of finite size. Nevertheless, as we discuss in the following, 
our technique still allows us to go beyond proof-of principle 
experiments and to simulate systems of a few tens of spins, 
which are no longer trackable using classical simulations. 

To proceed let us now for simplicity consider a homoge- 
nous ID system, which could be realized, e.g. using a tight 
ring trap. This system has periodic boundary conditions, and 
a constant lattice spacing a. The phonon spectrum can then be 
obtained in a closed form [34] and be written as 



tl<jj„ 




(18) 



where to a good approximation a>„ 2 Vl2| sta(m/N)\ is the 
dimensionless phonon frequency and y = Dm j (fid) » 1 is 
the ratio between potential and kinetic energy. Note that in 
contrast to the harmonic trap, here we assume mode labels n 
from —N/2 to N /2 and modes with identical \n\ are degen- 
erate. The maximum frequency of the phonon spectrum is 
the Debye frequency which in units of the dipole-dipole in- 
teraction D/(ha 3 ) is given by Cjd = u>„=n/2 — 6.94/ yfy with 
u>„ = <b n j y/y. The lowest frequency corresponding to the ro- 
tation of the whole ring is zero and doesn't play a role in the 
following discussion. We can write the normalized position 
fluctuations Sxj = Sxj/a as 



Oxjn/N + a \ e -i2njnjN 



), (19) 



from which we obtain a normalized spin-phonon-coupling 



8n,j = -i ■ 



-e 



ilnjnjN ~ 



2Nyaj„ 



(20) 



where g„ 6 sin(27rn/A0- Then the bare coupling = 
l/(2|z - j\ 3 ) and by defining w = yfyLL>/(D/ha 3 ) the phonon- 
mediated spin-spin coupling can be written as 



n>0 



gjcos[2jT(i-j)n/N] 
N(a> 2 - o%) 



(21) 



Note that the resulting phonon mediated couplings are inde- 
pendent of y, but a value of y » 1 has to be assumed for a 
(quasi) crystalline phase and the validity of our model. 

In Fig.|6]we plot the resulting total spin-spin couplings G,y 
for the case of N = 21 and different choices of the modulation 
frequency a>. For a low frequency a> = (a>i + a>2)/2 we obtain 
slowly varying, long-range interaction. Since the coupling of 
each spin to several nearby spins is positive and of similar 
strength, spin frustration effects as described above will also 
lead in this configuration to highly entangled grounds states. 
This is also true for the second example at a> — (a>3 + 0*4)12 
where the coupling is oscillating more rapidly, but still ex- 
hibits the properties G, , + i ~ G,-.,+2 > as well as £ &ij > 0, 



a cj — (a>i + tt>2)/2 
■ lj = (a>3 + U>4)/2 




8 9 10 



FIG. 6. (Color online) The normalized spin-spin coupling strength 
Gjj for N = 21 polar molecules in a ring trap as a function of the 
distance when we fix the dipole moment frequency at (a)\ + U)2)I2, 
(a>3 + £J<t)/2 and (a> g + a> i0 )/2, respectively. 



which prevents a simple anti-ferromagnetic or ferromagnetic 
ordering. Finally, the last example to = (<dc) + a>io)/2 shows a 
case where we address a high frequency phonon mode. Here 
the coupling exhibits site-to-site oscillations and we can ex- 
pect an anti-ferromagnetic alignment of spins in the ground 
state. Once the modulation frequency is tuned above a>o, 
the couplings start to decay faster until only the bare cou- 
pling G/j * G-9-/2 is left for detunings much larger than the 
mode spacing. Note that while in these three examples the 
total coupling to other spins 2 G,y is positive, also to oppo- 
site regime could be achieved by slightly detuning u> closer to 
one of the phonon modes. This would allow to tune between 
ferro- and anti-ferromagnetic long-range couplings. 

We now estimate the maximal number of spins that we 
could control by our technique. The validity of our model 
requires first of all that the effective spin dynamics is slow 
compared to the minimal detuning, i.e. \Gij/H\ < min{|A„ = 
cj n -a>\}. This can in principle always be achieved by choosing 
a lower value for e, but restricts the achievable spin-spin in- 
teractions to about ~ D/(Na 3 ). For a nearest neighbor dipole- 
dipole interaction strength of about 50 kHz and typical co- 
herence times in optical lattice experiments in the range of 
100 ms, this would still allow the simulation spin systems of 
N ~ 100 spins. 

A second condition for our model to be valid is that the 
maximal relative displacement Ax = max{|;t:,' - Jt, + i|} between 
the molecules remains small compared to the lattice spacing 
a. To evaluate if this condition can be satisfied, it is conve- 
nient to switch to a frame rotating with the modulation fre- 
quency a> and change into an interaction picture with respect 
to Yjn fiCL>ala„. Then the relevant slowly varying part of Eq. (0) 
becomes 

H = Zj hA " a n a " + Zj 2 8n ' ia " + 8n " 1 + ^ _ 2 -<T ' cr ''' 

n n,i i< j 

(22) 

In this form the Hamiltonian can be diagonalize by the 
unitary displacement operator U = exp(2„ o~(g* n ( .a„ - 
gnjcin) /2A„) 115 311 which displaces the phonon modes by a„ — > 
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a n - TjjSn ; <T z/(2A„), and correspondingly shifts the position 
of each molecule by 



N/2 



dXi 



2n{j - i)n 
N 



en 



(23) 



Since the system is homogenous we can chose i = and eval- 
uate the relative displacement \Sxq - Sxi\. Assuming that the 
main contributions come from the lower part of the phonon 
spectrum we can write the result approximately as 



N/2 



Ze g„ 2nn \ 
<■/(.. A N I 



w„A„ N 



(24) 



where5„ = ^ 2/ cos (~at) - 1- The relative displace- 
ments depend on the spin configuration S„. For a purely 
ferro- or anti-ferromagnetic ordering only S n= o or S„ = n/2 are 
non-zero, but for both modes g„ = 0. If we assume a com- 
pletely random distribution of spin states we can use the 
bound \S n \ < 1/ V^V and as above we consider a modulation 
frequency a> = (u) no + a> no+ i)/2 between two phonon modes. 
Then, apart from a numerical constant 0(1) we obtain the scal- 
ing l^io - Sx\\ ~ e V^V. Therefore, compared to the adiabatic 
condition mentioned above the assumption of small relative 
displacements does not impose any additional restrictions. 



VI. MOLECULAR SPIN QUBITS 

We finally come back to one of the initial assumptions made 
in our derivation and discuss in more detail, how to choose an 
appropriate encoding of the spin states \g) and \e) to obtain the 
correct form of the induced dipole moments yU gjf ,(f) as shown 
in Fig.[TJ(b). In the following we illustrate a basic encoding 
scheme for the case of a molecule with a 2 X ground state 15411 . 
where states \g) and \e) are represented in two Zeeman sub- 
levels of two excited rotational states. However, we point 
out that recently several other schemes for tunable dipole- 
dipole interactions for spin and/or rotational degrees of free- 
dom i37l l42l l43tl have been suggested which could also be 
adapted for the present purpose. 

We consider a 2 2 molecule in the vibrational ground state 
with a single unpaired electron with spin S — 1/2 and no 
nuclear spin. In the presence of a external field E = Ed c e z the 
rotational spectrum is described by the Hamiltonian 



H 



M 



BJ 2 -// c -E + y r J-S. 



(25) 



The first term is the rigid rotor Hamiltonian with B the rota- 
tional constant which is typically in the order of several GHz, 
and J is the total angular momentum. We denote the energy 
eigenstates of this rotor Hamiltonian by \J, M) with the spec- 
trum Ej = BJ{J+1), where M is the quantum number of asso- 
ciated with J z . The second term is the Stark interaction which 
occurs when an electric field is applied to a molecule possess- 
ing a permanent electric dipole moment // c . The external field 
mixes different rotational states and splits the (27 + l)-fold 



-0.08 

-0.10 

u -0.12 
a. 

4-0.14 

^-0.16 

v 

-0.18 
-0.20 







> 




^0 





2.0 2.5 3.0 3.5 

EJ(B/(i c ) 



4.0 



FIG. 7. (Color online) Induced dipole moments for |l,0)£ dc 

(dashed line) and \2,0)E ic (solid line) as a function of the applied 
dc field £d c . At the point £ , ( j c /(S//j £ .) = 3.05, those two states have 
the same induced dipole moment yu . In the range of rectangle, the 
dipole moment is approximately linear with external field, and the 
derivative of the dipole moments of two spin-states with respect to 
external field is the same but with different signs. 



degeneracy in the rotor spectrum, which amounts to induc- 
ing a finite dipole moment in each rotational state. We denote 
the rotational eigenstates in the presence of the bias field as 
\J, M) Eic . The third term in Eq. d25l ) is the spin-rotation in- 
teraction with y IS is the spin-rotation coupling constant. Typi- 
cally y rs ~ 100 MHz «B. 

For moderate and strong electric fields » B//i c , 
the Stark interaction typically exceeds the spin-rotation cou- 
pling and eigenstates of Hm are approximately given by 
product states \J,M) EiSc ® \M S ) Q with M s = ±1/2. 
Therefore, the dipole moment of each eigenstate is {M$ \ ® 
E6c {J,M\»z\J,M) Eic 9 \M S ) - Ed AJ,M\fi z \J,M) Eic . We encode 
our effective spin system in the two states \g) = 1 1 , 0) E ^ ® 
1-1/2) and \e) = \2,0) Ed . ® |l/2>, so that /i e = (e\fx z \e) = 
£dt <2,0k,|2,0> £dc and^ = <g|^|*> = £dc <l,0|^|l,0) £dc . For 
the spin-forbidden transition, transition matrix element be- 
tween \g) and \e) are small IB4I1 and compared with the in- 
duced dipole moments, the transition dipole moments (e\fx z \g) 
are suppressed by the ratio [y rs /B) 2 . However, the dipole- 
dipole interaction can still induce a spin preserving flip-flop 
process where the state of first molecule changes as 1 1 , 0) ® 
|-l/2> -> |2,0) ® I - 1/2) and the state of the second 
molecule as |2,0) ® |l/2> |1,0> ® |l/2> and brings the 
molecules out of the initial two level subspace. These pro- 
cesses can be suppressed by energetically shifting e.g. the 
state |2,0)®|-/2) -> |1,0)®| 1/2) by coupling it off -resonantly 
with a circularly polarized microwave field to the 7 = 3 man- 
ifold. 

In Fig. [7] we plot the induced dipole moments for two rota- 
tional states [l,0)£ dc and |2, 0) Etk as a function of the external 
electric field £dc- At the "sweet spot" E- dc - Eq = 3.05(B/fi c ), 
the induced dipole moments for those two states are both 
equal to /iq = -0.16/i c . In the range indicated by the rectangle 
in Fig. [7] we can make two approximations about spin-states 
\g) and \e). i) Their dipole moments are approximately lin- 
ear with the electric field strength; ii) The derivative of the 
dipole moments of two spin-states with respect to external 
field is approximately the same, but with different signs, i.e., 
d E 6c H s \E ~ -d E ^nA Eo = d Edc /i\ Eo . These conditions imply 
that when we fix the dc bias field at the sweet spot and add 
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another time dependent weak field £' ac (f) = £ ac cos(Lot), the 
dipole moments of two spin states become fi{t) = yUo + Snifter 1 
with 6fi(t) = 3£ dc /i|£ £ac cos(w/)- Note that within the linear 
regime of the dipole moments variations of e — \Sfi\/fXo ~ 0.1 
as required for our model can be achieved. 



VII. CONCLUSION AND OUTLOOK 

In this work we have studied the interplay between direct 
and phonon-mediated spin-spin interactions in dipolar crys- 
tals of polar molecules. In particular we have shown that 
for an appropriate choice of rotational states an electric bias 
field with time varying amplitude can be used to induce a 
state dependent modulation of the dipole moments of the 
molecules, which can resonantly enhance the coupling of the 
internal degrees of freedom to a specific collective phonon 
mode in the crystal. This minimal level of control already 
enables the implementation of Ising models with non-trivial 
spin-spin interactions which for small and moderate sized sys- 
tems can be used to generate various multi-particle entangled 
states or to simulate Ising models with 'infinite-range' inter- 
actions. To extend this method to larger spin systems and 
to achieve a more flexible control over the resulting spin- 
spin interaction patterns different generalizations of our tech- 
nique can be envisioned. This includes arrays of microtraps 
where within each trap a small number of molecules inter- 
act via phonon-mediated couplings while couplings between 
traps are achieved via direct dipole-dipole interactions. In this 
case the frequency splitting between phonon-modes does not 
decrease with the total systems size. Further, in analogy to 
Ref. 114311 , local defects or so-called 'marker' molecules can 
be introduced to modify the local phonon structure and by that 
tune between long- and short-range interactions. 



Appendix A: Effective spin dynamics and state purity 

In the absence of any transverse magnetic fields the full 
time dependent Hamiltonian given in Eq. (0 can be integrated 
exactly and the resulting time evolution operator can be writ- 
ten as 



Here the displacement amplitudes are 



a n ,i = i- 



l — e -<(0Jn-0>)t I _ 



U> + Cl>„ 



(Al) 



(A2) 



and the phases 0, ; (f) = <X>y(0 + $>h(t) are given by the bare 
coupling 



Go 
(1 a sin(2wf) 

and the phonon-mediated part 

2g n ,ign,j(*>n r t sin(2wf) 



(A3) 



u> n cos(wf) sin(w„f) + co sin(wf) cos(w„f) i 
{co> -aft J 



(A4) 



From these expressions we can extract the slowly varying 
spin-spin couplings G-'; 1 := lim,^ Of; 1 /t. 

To evaluate the purity of the spin subsystem during the evo- 
lution we write an arbitrary initial spin superposition state as 
l^o) = Zs>0|s) where cr' z \s) = Si\s). With p ph being the equi- 
librium density operator of the phonon modes the total density 
operator at time t — is p(0) = lYoXYol ®Pph and the reduced 
spin density operator at time t is given by 



p s {f)=Yu c ^y m,lf)e ' 1 



-kZnFnitA?)] 



MA- 



(A5) 
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Here 0(f, s,r) - Xi<j ~ nrj) and 



F n (t, s, r) 



^Onjim-?,) coth|^j, (A6) 



where T is the temperature of the phonon modes. The purity 
P s (t) = Tr{p 2 (f)} of the spin system can then be written as 



p,(0 = Xiqi 2 ic;iV2»^)^. 



(A7) 
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